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1 Introduction and Summary 



The Maldacena Conjecture p] provides what are probably the most effective 
and controlable tools to study non-perturbative dynamics of a variety of field 
theories. A large variety of effects have been discovered or checked in field 
theories using a suitable string dual. Integrability, correlation functions of 
various interesting operators (protected or not by symmetries), aspects of 
lower dimensional systems, applications in condensed matter and QCD-like 
systems have been succesfully studied using gauge/gravity duality. 

The results above, while more numerous and spectacular in highly (su- 
per)symmetric theories, are not restricted to examples of this sort. As a 
matter of fact, there are many applications where black holes (and hence 
dual field theories at finite temperature) play a fundamental role. In these 
cases the dynamics is neither driven by SUSY nor by conformal symmetry. 

As a result, an interesting problem is to construct backgrounds duals to 
field theories where supersymmetry has been broken in a soft way. These 
systems should conserve some of the dynamics of the SUSY case with the 
addition of the deformations by relevant operators that break the supersym- 
metry. The low energy dynamics should then be a non-linear superposition 
of SUSY and non-SUSY effects. This is an interesting problem, on which it 
seems feasible to make progress. 

In this paper, we will construct duals to field theories in four dimensions 
where SUSY has been explicitly and softly broken by the addition of rele- 
vant operators to the Lagrangian. The original field theories will be those 
obtained by a twisted compactification of five branes wrapping a calibrated 
two cycle in the resolved conifold and those obtained by studying the dy- 
namics of D3 and fractional D5 branes on the tip of a conifold. Both are 
non-conformal theories with interesting low energy dynamics (confinement, 
R-symmetry breaking, formation of domain walls, k-strings, etc.) 

We will construct our non-SUSY backgrounds by finding an explicit so- 
lution of the Einstein, dilaton and RR-form equations of motion. We also 
impose that irrelevant operators are absent from the dynamics and that the 
string backgrounds are regular all along the space. We will concentrate on 
the case in which the SUSY breaking parameters are small compared the 
others already present in the system in the SUSY case. 
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These will then be examples of backgrounds dual to the strongly coupled 
dynamics of well understood SUSY field theories in which SUSY has been 
softly and controllably broken. Some examples of this sort have appeared in 
the past for deformations of well-known SUSY backgrounds, see for example 

Our paper is organized as follows. We start in Section [2] by presenting 
the SUSY system. While the fomalism summarized there does not apply to 
problem of interest, we do give some details that are useful in attempting to 
construct the non-SUSY solutions (in particular large-radius asymptotics) . 
In Section [3j we will the propose a SUSY breaking solution in series ex- 
pansion for large (UV) and small (IR) values of the radial coordinate. We 
will carefully count the parameters that control our solutions and find nu- 
merical solutions interpolating in a smooth way between the desired UV 
and IR asymptotics. Section |4] gives some details of the numerical method. 
In Section [5] we calculate the ADM Energy of the new solutions (with the 
SUSY solutions as reference backgrounds). In Section |6] we perform a de- 
tailed study of various field theory quantities, whose strong-coupling result 
points us to an interpretation of the dual field theory being deformed by 
the insertion of relevant operators, like gaugino masses that break SUSY 
and may also infiuence VEVs. We close with some conclusions and possible 
interesting problems to be solved constructing on the results of this paper. 
The high technical nature of our work is clear from the outline above. For 
the benefit of readers, we have included explicit technical points in detailed 
appendices. 

Note Added: While this paper was close to completion, we were informed 
of the work by Dymarsky and Kuperstein, having interesting overlap with 
ours ]35^ . We thank Anatoly Dymarsky for letting us know about this work 
prior to publication and discussion on these topics. 

2 Presentation of the SUSY system 

In this section we summarize well established aspects of particular super- 
symmetric field theories and their dual backgrounds. This will be useful 
when introducing SUSY breaking deformations. 

We start by considering two apparently different field theories. The first 
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one, we refer to it as 'theory I' or 'Type I theory' (hoping not to cause 
confusion with the Type I string theory), is a quiver with gauge group 
SU{n+Nc) X SU{n) and bifundamental matter multiplets Ai^ Ba with i,a = 

I, 2. The global symmetries ar^ 

SU{2)l X SU{2)r X U{1)b x U{1)r. (2.1) 
These bifundamentals transform under the local and global symmetries as 
Ai = {n + Nc,n,2,l,l,^), = (n + iV„ n, 1, 2, -1, ^). (2.2) 
There is also a superpotential of the form W = -(-ij^ap 

tr[AiBaAjBp\. The 

field theory is taken to be close to a strongly coupled fixed point. In that 
case one can show that the anomalous dimensions should be 7a,b ~ —\- 
This field theory is well known to be the dual to the Klebanov-Strassler 
background [3] and its generalization to the baryonic branch [3]. 

The second field theory, that we will call 'theory IT (again not to be con- 
fused with the Type II string!) is obtained after a twisted compactification 
(to four dimensions) of six dimensional SUSY SU{Nc) Yang-Mills with 16 
supercharges. This special compactification studied in [5], [6] preserves four 
supercharges. In four dimensional language, the field content is a massless 
vector multiplet and a 'Kaluza-Klein' tower of massive chiral and massive 
vector multiplets. The Lagrangian, the weakly coupled mass spectrum and 
degeneracies are written in [6]. The local and global symmetries are (the 
R-symmetry is anomalous, like in the theory I above), 

SU{N,) X SU{2)l X SUi2)R x U{1)r. (2.3) 

These two theories, apparently so different, can be connected as discussed in 
[7] and [H] via higgsing. Indeed, giving a particular (classical) baryonic VEV 
to the fields {Ai, Ba) and expanding around it, the field content and degen- 
eracies of [6] is reproduced. This weakly coupled field theory connection has 
its counterpart in the type IIB solutions dual to each of the field theories. 
Indeed, it is possible to connect the dual backgrounds to field theories I and 

II, using U-duality [7J. This connection was further studied in [8], f9], [TOj, 

m- 

^The R-symmetry is anomalous, breaking U{1)r Z^n^- 
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We will now explain this connection among the explicit Type IIB string 
backgrounds. We start from a background describing the strong dynamics 
of the 'field theory II' (the twisted compactification of five branes). A quite 
generic configuration of this kind can be compactly written using the SU{2) 
left-invariant one-forms 

ioi = cos ipdO + sin ip sin 6d(p , 0)2 = — sin Tpd9 + cos ^p sin 9d(p 
= dip + cos9d(p (2-4) 

and the vielbeins 

= e^dxi, EP = e^+''dp, E'^ = e^+^d9, E'^ = e^^^ sin9d^, 
^1 ^ ie?+»(wi + adO), E^ = ^e^+^C^a - asm9d^), 

= ^e?+''(LD3 + cos ed^). (2.5) 

In terms of these, the background and the RR three-form read 
10 



dsl = Y,{E^)\ (2.6) 

i=l 

F3 = [/i^^^^ + /2^^^^ + HE^"^^ + E^^^) + fi{EP^'^ + EP^'^) 
where we defined 

Eijk..l A E^ AE'' A ... A E\ 

h = -2N,e'^-^\ h = ^e-'~'\a^ - 2ab + 1), 

f, = N,e-^-^-^{a-b), U = ^e-^-^~%'. (2.7) 

The dilaton, as usual, is a function of the radial coordinate ^{p) and we 
have set a' Qs = 1. 

The full background is then determined by solving the equations of mo- 
tion for the functions (a, 6, 5, h,k). A system of BPS equations is derived 
using this Ansatz (see appendix of reference [12] ) . These non- linear and cou- 
pled first order equations can be arranged in a convenient form, by rewrit- 
ing the functions of the background in terms of a new basis of functions 
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P{p), Q{p),Y{p), t{p), a{p) that decouples the equations (as explained in 
[13] -[H]). We quote this change of basis in our Appendix [A| 

Using these new variables, one can manipulate the decoupled BPS equa- 
tions, solving most of them and obtaining a single decoupled second order 
equation for P{p)- All other functions are obtained from P{p) — see [13] and 
our Appendix [A] for details. The second order equation mentioned above 
reads 



We will refer to Eq. (2.8) as the master equation: this is the only equation 
that needs solving in order to generate the large classes of solutions of Type 
IIB dual to "field theory 11" in different circumstances (vacua, insertion of 
operators in the Lagrangian, etc.) 

In this paper, we will not be concerned with SUSY solutions, but they 
will play an important guiding role. We summarize below the small and 
large p expansions of the function P{p). 

2.1 Aspects of the SUSY solutions 



Let us start from the solution of the master equation (2.8) for large values of 
the radial coordinate (describing the UV of the field theory II) . The SUSY 
solutions have an expansion for p — t- oo of the form. 



P = e^Pl^ 
^4g-i6p/3 / 18567 2781 27 



„2 QR„3 

c\ V 512 32 ' 4 ' 



+ 1 + -tP - 36/0' 



(2.9) 



Notice that this expansion involves two integration constants, c+ > and 
c_. The background functions at large p are written in Appendix [Xj 

Regarding the IR expansion, we look for solutions with P — )• as p — )• 0, 
in which case we find 



4/ii / 4iV2\ ig;^ / 32Ar4\ 



(2.10) 



^As an example, the solution P — 2NcP gives the background of [5], [IS]- This solution 
and those with the same large p asymptotics will not be the focus of this paper. 
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where hi is again an arbitrary constant, there is of course another integration 
constant, taken to zero here, to avoid singularities. This gives background 
functions that are quoted in Appendix [Aj Of course, there is a smooth nu- 
merical interpolation between both expansions. However, there is then only 
one independent parameter; given a value for one of {c+, c_, hi}, the require- 
ment that the solution matches both expansions is sufficient to determine 
the values of the other two. 

As explained in [SJ, this solution corresponds to a dual field theory II in 
the presence of a dimension-eight operator inserted in the Lagrangian which 
ultimately couples the field theory to gravity. This calls for a completion in 
the context of field theory. This is achieved with the U-duality of [7] (which 
we will sometimes refer to as the 'rotation'). 

After the U-duality described in [7] is applied, we define the new vielbein 
(which we use in the following). 





e^h 


e' = 


1 * 




1 # 
2^" 



* 

64' 



'hidp, 



64 



6T+^^/i4(tii +ad0), 6^ 
'^''h^ojs + cos 9dip). 



1 $1^1 

-e^ ^hi {uo2 — a sin 9dip), 



hi sin Odip, 



(2.11) 



The newly generated metric, RR and NS fields are 

10 

dsi = ^(6^)^ 



i=l 



[/ie^23 + f2e 



e<^3 ^ ^.^(^623 ^ g^l3) ^ y^^gpie ^ ^p^2. 



64 



-K 



/j3/4 

K^^dt A dxi A dx2 A dx^, 
h 

2*^ 



193 _|_ gi^23> 



= Ke 4 
We have defined 



'h-id. 



h 



(2.12) 



h = l- K'e"", (2.13) 

where k is a constant that we will choose to be k = e~^^°°\ forcing the 
dilaton to be bounded at large distances. The rationale for this choice is to 
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obtain a dual QFT decoupled from gravity. Details of this were carefully 
discussed in [8], [11]. The tuning k = e~*(°°-' (also chosen in [16j, though 
in slightly different notation) is the geometric version of the fact that, in 
order to eliminate an irrelevant operator in the dual field theory I, we have 
to finely-tune the matter content and the gauge group with which we will 
UV-complete the theory II after un-higgsing from the single node to the 
quiver. See [8] for a complete explanation. We will now move to study 
SUSY breaking deformations. 



3 The SUSY-breaking deformation 

The goal is to find a non supersymmetric solution with the same symme- 
tries and structure as the ones described above. We proceed as follows: 



we will find a non-SUSY generalization of the system in eq.(2.6). We will 
solve the equations corresponding to Einstein, Maxwell, dilaton and Bianchi 
equations of the system. The nice properties of the SUSY formalism just 
explained do not apply. We will then propose series expansions for the in- 
dividual background functions h, g, k, o, b. With the experience gained in 
the SUSY example, especially keeping in mind the expansions quoted in eqs. 



(A. 3 A. 5), we propose similar asymptotics. 
3.1 Asymptotic expansions 

In the UV (large values of p) our expansions take the form, 

OO 2 2g ^ 



-i)p/3 
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i=0 j=0 1=0 j=0 

2/j. OO 2 OO 2 

i=0 j=0 1=1 j=0 

OO z OO i 

bip) -EE^^^f^^"^'-'^'^"' <p) -EE^^^p^^""-""^'"- (3-1) 

i=l j=0 i=l j=0 

We have found that a generic solution of this sort can be written in terms 
of nine integration constants. These constants are free; all other coefficients 
in the series expansion can be written in terms of them. The independent 
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constants are taken to be, 

Kqq, K30, Hio, Hu, ^10, ^30, W20, W40, V40. (3.2) 

Note that we have found the constants V21, W21 must vanish for this to be a 
solution. Also, we imposed that terms that would spoil the UV behavior of 
the SUSY solution (corresponding to irrelevant operators in the dual QFT) 
are absent from our expansions. 

Without loss of generality, we relabel the UV parameters in eq.(3.2) to 
make contact with the SUSY case (see appendix |A]) : 

2c, 



,3 



The independent parameters are then 

c+, c_, $00, Qo, Po, Hu, W20, ^30, V40, (3.4) 
and we can recover the SUSY case by setting 

^11 = 7T, W20 = 0, $30 = 7^(3 + 4Qo), F40 = 2e2''°(l + Q„). (3.5) 

For small values of the radial coordinate (which we take to end at p = 0) 
we will propose an expansion of the form (again, imposing regularity of the 
solution), 

e''~>>.y, ^-r-y:9,f^, ^~>>.y, (3.6) 



4$ 

e ~ 



In this case the free parameters are /cq, /o, k2, V2 and W2- For any value of 
these numbers we find a solution. To make contact with the SUSY solution 



00 


^ 


g2fc °° 

— 




j=0 


j=0 


00 


00 

o(p) ^^Wjp>, 
j=0 


00 

Kp) ~ E ^^p' 

j=0 



in eq.(A.5), we relabel = hi/2 and /o = e^"^". We then recover the SUSY 
solution if the remaining three parameters take the values 

,, = _L(,;_4), = a., = ±-2. (3,7) 
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If we want to restrict our attention to those solutions which match both 
the UV and IR expansions described above, we expect to have fewer inde- 
pendent parameters. To see this, note that we can describe the solutions 
either by the IR boundary conditions, so that they are parameterised by the 
five IR parameters {hi,4'Q,k2,V2,W2}, or by the UV boundary conditions, 
giving a parameterisation in terms of 

{c+ , C_ , $00 , Qo, /Oo, -ffii , W20, $30 , V40} • 

If a solution exists which connects our IR and UV expansions, the func- 
tions resulting from these two parameterizations must be the same. We can 
formally express this as a system of twelve equation^ 

ghi...w2{p) = 9c+...vM, 9'hi...w2(p) = 5c+...i/4o(^)' 
K^...W2{P) = hc+...vM^ K^...wM = K+...Vao^P)^ 

K-^2ip) = bc+...vM^ K^...w2ip) = ^'c+...v,M- 

However, the derivative of one of the functions can be expressed in terms 
of the other derivatives and the functions themselves using the constraint, 
so only eleven of these equations can be independent. We therefore expect 
to be able to solve for eleven of the fourteen parameters, leaving only three 
independent. 

We know that the dilaton can be shifted without otherwise affecting the 
solution, so one of these parameters must be either (j)Q or ^oo- Additionally, 
we know that we also need one of hi, c+ or c_ to describe the class of SUSY 
solutions. The third parameter therefore breaks SUSY. 

Note that there does not appear to be anything to stop us choosing, say, 
Qo to parameterise the SUSY-breaking, despite the fact that UV expansions 
with Qo —Nc do not break SUSY. This is explained by the fact that SUSY 
solutions with Qo 7^ —Nc do not have a regular IR. If we simultaneously 



demand that Qo 7^ —Nc and the IR is of the form (C.13-C.24) we must 



therefore have a non-SUSY solution. However, it seems conceptually simpler 



^We write the functions resulting from a given choice {/ii , 00, ^2, 1^2, ^"2} of the 
IR parameters in the form 5'hi,0o,fe2,f2,ii'2 (p)- Similarly, the expressions of the form 
5c+,c_ ,*oo,<3o,pe.,irii,W2o.*3o>V4o (p) Tsfer to the functions resulting from a given choice of 
the UV parameters. 
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to choose the third parameter to be one which explicitly breaks SUSY. For 
our three independent parameters, we might select (in the IR) {hi,(j)o,W2}, 
or (in the UV) {c+, $00, ^20}- 

In the next section, we wih study the challenging numerical problem 
of finding a solution that interpolates between the small and the large p 



expansions in eqs.(3.1 )-(3.7) 



To summarize: we want to find a numerical solution for the functions 
above. This will provide us with a non-SUSY deformation of the background 



in eq.(2.6) dual to field theory of type II. Applying the U-duality in [7]-[TT] 
we construct a non-SUSY background of the form given in eq.( 2.12[ ), dual 
to a non-SUSY version of the field theory I. 

4 Numerical analysis 

In this section we show that for some values of the free constants in the IR 
and in the UV we can connect the asymptotics numerically. We first briefiy 
describe our method of relating the IR and UV parameters, the details of 
which are relegated to appendix|D| before presenting a sample solution which 
smoothly connects the IR and UV asymptotics of the previous section. 



Our approach is to solve the equations of motion (B.4-B.9) numerically, 
using the expansions (see appendix [C]) as boundary conditions. We start 
from the IR expansions ( |C.13 -C.24), meaning that our numerical solutions 



are described by the three SUSY-breaking parameters {k2,V2,W2}, in addi- 
tion to those already present in the SUSY case hi and (pQ. 



However, we have seen in section 3.1 that we expect only three inde^ 



pendent parameters in total, one of which breaks SUSY. This would mean 
that in the non-SUSY case we cannot treat even the five IR parameters 
as independent. This is confirmed by our numerical analysis: a generic 
choice of the IR parameters yields UV behaviour of the form b ~ ite^^ and 
g29 g2/i ^ g2fc ^ e^p/^^ which is incompatible with our expansions (3.1). 



To obtain a solution connecting the IR and UV expansions, then, we have to 
determine both an appropriate combination of values for the IR parameters 
and the corresponding values for the UV parameters. 

We achieve this using the method described in detail in appendix |D] In 
outline, we start with a manual search of the IR parameter space, using 
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Mathematica's NDSolve to obtain numerical solutions. Having obtained a 
solution with the desired UV behaviour (b 

we optimise the match to the UV expansions (C.1-C.6) with respect to the 
parameters (3.4) using NMinimize. 

An example of solution is shown in figure [TJ As is expected from the 
expansions (C.IC.G), the most significant modification with respect to the 
SUSY solution is the presence of the e"^''/^ behaviour in the UV of a and b. 
The size of this effect is controlled by the SUSY-breaking parameter W20', in 
the SUSY case we have W20 = 0, resulting in a ~ e"^'' and b ~ pe^'^^. The 
other functions are modified at higher order in (C.l C.6), and as expected 
no effect is visible in figure [T] See Appendix [D] for details of the numerical 
analysis. 



5 Energy 

In this section we study the energy of the non-SUSY solutions found above. 
For any stationary spacetime admitting foliations by a spacelike hypersur- 
face Si, the free energy and the energy are related via the thermodynamic 
relation F = E — TS. Here we are considering T = backgrounds and so 
we expect F = E. In this section we will first calculate the ADM energy 
E, for the solutions before the U-duality — we will refer to the U-duality 
of reference [7] as 'rotation'. We will then repeat this calculation for the 
solutions after rotation and show that the energies before and after rotation 
are equal. As a check of our results, in Appendix [E] we obtain the free energy 
using the on-shell action method and show that F = E. 

5.1 ADM energy 

Consider a non-asymptotically flat 10-dimensional background. Let T,t be a 
9-dimensional constant-time slice whose 8 dimensional boundary is a constant- 
radius surface Sf°. The regularized internal energy E is defined as [T7] . 




Nf is the lapse function, Nj^ is the shift vector, p^jy the momentum conjugate 
to the time derivative in the constant time-slice, and ^Kq are the extrinsic 
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p 



p 



5 10 15 5 10 15 

log a log b 




P 



Figure 1: Plots of the functions (7, /i, fc, log a and log 6, obtained numeri- 
cally (solid blue), together with the IR (dotted red) and UV (dashed orange) 
expansions (appendix [C]), with small deviations from the SUSY values of the 
parameters. The SUSY solution (grey) is included for comparison. 
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curvatures of the 8 dimensional boundary S^, for the background under 
consideration and the reference background respectively. Finally n'^ is the 
spatial unit vector normal to the constant radius-surface Sf^. It is required 
that both geometries induce the same metric on St°°. The matter fields 
should also agree at Sf° or at least the difference should tend to zero as 
goes to infinity. We will choose a SUSY background as a reference geometry. 



$/4g The extrinsic curva- 



For the metrics before rotation ( 2.5 - 2.6 ) we have Nj^ = 0, Nt = y^j^oo 
ture is 

1 



-$/4-fc 



2{^' + g' + h') + k' , (5.2 



where gg denotes the determinant of the 9-dimensional constant time slice 
Sj. The requirement that the induced metrics on St°° agree at the boundary 
impliean 



e 2 e 



23n 



e 2 e 



29s 



e 2 e 



2h„ 



e 2 e 



2hs 



e 2 e 



e 2 e " 



and the goo component agrees if 



e 2 = e 2 . 



(5.3) 



(5.4) 



All the quantities in (5.3) and (5.4) are evaluated at some large but finite 



rc that acts cutoff. Using (^5.1|, the energy is 
vols 



E 



647r 



lilll ^ 6 ^iT'S ( -|-2^7ts -|-2/i7is -|-A;t) 



g-fcs / ^2<S>s+2gs+2hs+ks 



(5.5) 



Before evaluating (5.5) we have to satisfy the matching conditions at the 



boundary, (5.3) and (5.4). In order to do this we have to use the most 



general asymptotics of a supersymmetric solution. As discussed in eq.(3.5), 



analyzing the BPS equations (see Appendix A in reference [12]) we see that 
the most general supersymmetric UV asymptotics is obtained by replacing 



20 



0, ^40 ^ 2e2^°(l + Qo), Hu ^ 1/2, $ 



30 



(3 + 4Qo) 



4d 



(5.6) 



*The subscripts ns and su stand for non-supersymmetric and supersymetric respec- 
tively. 
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in the non-supersymmetric expansion (3.1). Notice that this substitution 
restores the integration constants Qo, Po and e*°° that are usually set to 
— 1, and 1 respectively [13J. Reintroducing the integration constants is 
equivalent to using the shift invariance of the r coordinate (encoded in Qo 
and Po ) and the dilaton [18]. Adjusting these constants will allow us to 
satisfy the matching conditions at the boundary and cancel divergences in 
the energy. Given the complexity of the UV expansions the matching pro- 
cedure is cumbersome but straightforward. Working to linear order in W20 
we obtain, 

S= -^4e2^«+2*-VF20. (5.7) 
After the duality transformations the UV asymptotics changes drasti- 
cally. In this case we have =0, Nt = y Iffool = ^jpjT^ dSf° = 

lg3$+23+2/^+fc^l/2^ ^ g-3<I./4g-fe^-l/4 



K = V^n'^ = — 9,.(V9^n'^) = [H' + 2H{3<!>' + 2g' + 2h' + k')] . 

(5.8) 

Note that here we defined H = e~^* — e"^*°°. The regularized energy after 
the rotation is 

E = -^lim{Ar,,-Asu} (5.9) 

047r r— >oo 

where 



The matching conditions now read, 

Trl/2 ^+2g„s ^ rrl/2 ^+2gsu ttI/2 ^+2hns ^ tt1/2 ^+2hsu 

ifV2e^+2fc„. ^ ^^^^3^+2..„^ Hje-"^ = H^^'/h-'^. (5.11) 
Note that 

, s/ p<S>+2g+2h . _ • 



'H 

= /^before ^extra ^^ -^2) 

where A^'^/"- = ^e'^+29+2h+ky ^extra = -^^^ (c'^Vh)' . We 
have 
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where all the functions are evaluated at some large but finite cutoff rc- 
After adjusting the parameters to ensure that the induced metrics at the 
boundaries are the same, as required in (5.11 ), we take the cutoff to infinity. 
The first two terms in (5.13) are the same as in the energy before rotation 
(5.5). We find that — to first order in W20 — the matching conditions are 
satisfied using the same set of integration constants as before the rotation. 
Thus, the first two terms in (5.13) give exactly the energy before rotation. 
Any difference in energies will come from the extra terms (5.13). However, 
it can be shown that using the integration constants necessary to satisfy 

( [5lT| ), 

lim {(A^f ™ - *™)| = 0. (5.14) 

Thus the energy before and after rotation are the sam^ Indeed, plugging 
in the UV expansions directly in (5.9) we obtain, 

-^4e2^«+2*-TF20. (5.15) 

A couple of comments are in order. First, note that the overall constant 
that appears in the energy can be changed by shifting the value of the 
dilaton at infinity. Thus, the physically meaningful statement is that the 
energies before and after rotation have the same functional dependence on 
the parameters, 

Second, this calculation can be carried out to higher order in the SUSY 
breaking parameter W2o- The divergences in the energy can be cancelled 
by subtracting an appropiate SUSY background. However, at higher orders 
there will always be a discrepancy of order VFfo of the metrics at the bound- 
ary. This clearly indicates that the treatment presented in this section is 
valid only for soft supersymmetry breaking with small breaking parameter, 
W20- Had we not expanded around W20 ~ the mismatch at the boundary 
could be arbitrarily large indicating that the non-supersymmetric solution 
does not approach the SUSY solution fast enough for the energy to be finite. 
Note that this substantiates the smallness of W20 seen numerically in the 
previous section where the solutions found have W20 ~ 0(10^^). 



^This suggests that the ADM Energy is 'uncharged' under the U-duaHty, hke probably 
are also uncharged various thermodynamical quantities. 
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6 Field Theory Aspects 



In this section we will analyze various field theory aspects of a non-SUSY 



version of the quiver that we called field theory I and described below eq.(2.1 ) 



To this end, we will use the non-SUSY background one obtains when plug- 
ging our numerical solutions in Section[4]in the background of eq.(2.12 ) dual 
to the field theory I. 



To begin with, notice that in eq.(2.12) we did not specify the NS potential 



i?2- Since this will be useful below, we discuss it here (the result is different 
from the SUSY one). 

Following the intuition gained in the SUSY example, we propose a B2 
of the form 



B2 = bi{p)eP' + b2{p)e'^ + bsip)e'^ + b^ip)e'^ + b5ip)e^\ 



(6.1) 



by imposing that dB2 = H3 and that the Page charge vanishes Qpage, D3 = 
(see below) we obtain — all details are discussed in Appendix |F] — 



b2 



4h 

,-2h 



263$' - 3hb3^' - Ahbsg' - 2hb'.^ + kN^c'^-'^^U [a^ - 2ab + l) 



4/ii/2 
bi = b^ = - 



e^9hl (1 - a") 63 



K 3* 

— e 2 

Nr. 



iV2(a- 6)6' + 46^(3+'^)$' 



1 . KN^e^~9-H' 
-^^"^063 



(6.2) 



2 4/ii/2 

with b-i{p) an undetermined function. This freedom corresponds to a gauge 
transformation. A general B2 can be expressed as 

B2 = (i?2).3=o - Id^e'^-'^^'^h^'^, e') . (6.3) 

Before computing various observables of the strongly coupled non-SUSY 
field theory I, we will quote another quantity that will appear frequently in 
the analysis. This is a periodic quantity in the string theory. Given the two 
cycle defined as, 

S2 = [0 = ^,99 = 27r-(^,V = V'o], (6.4) 

we define 



60 



1 



B2. 



(6.5) 
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When computed explicitly using the form of the B2 potential in eqs.(6.1 6.2), 
we obtain 

47r ttNc 



These quantities together with those appearing in the background of eq. ( 2 . 1 2 ) 



will be important in the study of the non-perturbative field theory dynamics. 
6.1 Calculation of observables 

We now move into the calculation of observables that will help us understand 
the field theory interpretation of our solution. 

6.1.1 Interesting Asymptotic Behaviors 

We start by studying some combination of fields that have a particular 
behavior. For this it is convenient to reduce the system to five dimensions 
as was done in j8]. Once in five dimensions, the paper [8] shows that some of 
the 5-d fields are invariants under the rotation. These fields are the dilaton 
^ and the combinations 

Ml = 1 + + 462^^-29, M2 = e2''+2s-4fc (5.7) 

The corresponding UV expansions are (see Section [4] and Appendix [C] for 
the notation) 

e-4p/3 

Ml = 2+ {8Hup + 3c+W^o + 2Qo) + 0{e-^P/^) , 

^+ 

= ^ - ^<e-'''/' + +0{e-^Pl^). (6.8) 

Now, suppose that we define a variable z = e"^''/^. Any field M that 
for z — 7- scales like ~ indicates either the insertion of a rele- 
vant / marginal operator or the VEV for an operator of dimension A (if A > 
or A = 0). On the other hand, if A < 0, it indicates the insertion in the 
Lagrangian of an irrelevant operator of dimension (4 — A). 

Using the UV expansion of the dilaton (see Appendix [C]) , we see that 
the dilaton corresponds to a marginal operator of dimension A = 4 (this is 
identified with a combination of gauge couplings g\ discussed below) . The 
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expansion of the function b{p) — see again Appendix [C| — indicates that the 
SUSY breaking constant W20 corresponds to the insertion of an operator of 
dimension three in the lagrangian. We associate this operator with the mass 
for the gaugino and in an analogous way, the constant e^''° is associated with 
the VEV for the gaugino. The association is not exact, in the sense that 
once SUSY is broken there could be a contribution of W20 to the gaugino 
VEV. Then, schematically we have 

W20 mAA, e'^P° (AA) ~ A^^- (6-9) 

Following this logic, the expansion of the field Mi ~ is interpreted as the 
VEV for a dimension two operator [16j, 

U ^ tr[AA^ - B]. (6.10) 

This same operator gets a VEV in the SUSY case and is the one that allows 
us to explore the baryonic branch. Notice that the SUSY breaking coefficient 
W20 contributes to this VEV. 

In the theory of type I, that has two gauge groups, we should expect two 
independent gaugino masses. Here, the solution is obtained by U-duality 
applied on a background dual to a Theory of type II, with only one gauge 
group. We seem to have only one integration constant associated with gaug- 



ino mass, that is W2o- As emphasized below eq.(3.2) the numbers V21, W21 
corresponding to behavior of the functions a ~ 6 ~ pe~'^P^^, which could 
be associated with this second mass parameter, turn out to vanish in our 
particular solution. 

6.1.2 Energy 

We take the expressions for the ADM Energy of the non-SUSY backgrounds 



as derived in eqs.(5.7) and (5.16) and we use the map described in eq.(6.9), 
we obtain that 

Eadm ~ 4e2*(-)e2''0T^20 ~ mAf.^. (6.11) 

Then the energy is proportional to the gaugino mass and the strong coupling 
scale, as expected. 
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6.1.3 Charges 

We will define the Maxwell and Page Charges 

^Maxwell, D3 = ^ n a / -^Sj ^Maxwell, D5 = / -^3) 



QPage, D3 = j ^5 - B2 A F3, 



(6.12) 



where the manifold X5 = [6, ip, 9, ip, ip] and X3 = [6, (p, Tp]. As in the SUSY 
case we have that 

^Maxwell, D3 = -e^^^^'^'^^*<I>', ^Maxwell, D5 = ^c- (6.13) 

vr 

We have also imposed that Qpage, Ds = in determining the B2 field of 
eq.(6.1 ) — see Appendix [F| for details. The vanishing of the D3-Page charge 



is a feature of the SUSY non-singular solutions; this is the reason why we 
imposed it here. It would be interesting to see if one can obtain a regular 
non-SUSY solution in the presence of sources indicated by a non-vanishing 
Page charge. Using the UV expansions, the Maxwell charge for D3 branes 
is 

e*- 1 



QMaxwell, D3 = " ^ (9e*°° + 4^6 ^*-$3o) 

33e*°°W"o^ 



327r 

So, we see that W20, the same number that determines the mass of the 
gaugino according the discussion above, changes the large energy value of 
the Maxwell charge (correspondingly of the c-function — see below) in a 
subleading way, as expected. 



6.1.4 Gauge couplings and beta functions 

Let us review briefly what happens in the SUSY case. In the SU{Nc + n) x 
SU{n) SUSY quiver, we have two couplings (71, 92- Close to the Klebanov- 
Witten conformal point (in the UV), the anomalous dimensions are ^a,b ~ 
— ^. This implies that the beta functions for the diagonal combinations 

/38.2 = /3«,2 - /38.2 = 6N„ (3,^2 = P,^2 + = 0. (6.15) 

9?r iif ^ ^ ^ 
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As in the SUSY case, we will adopt the definition^ 



^ = ^e-*, ^ = 2^e-* [1 - 6o(7r)] (6.16) 

9+ 9- 



where bo{'tpo) is defined in eq. (6.5)-(6.6). We obtain 
5- V / 2 



Notice that the result is independent of the gauge artifact function 63 (p) . In 
the UV, these formulas are typically trustable. The explicit expansions are 

(6.18) 

and 



_/(iT> „ _(T) <3 \ O 

(6.19) 

Let us now compute the beta functions as read from the geometry. We will 
use the radius/energy relation 

r = e^Pl^ = ^ (6.20) 

where is the energy scale at which we probe the process and A the ref- 
erence or strong coupling scale of the given gauge group. Notice that this 
choice is arbitrary, just refiecting the possibility of choosing a scheme. Other 
monotonic relations p(//) would express the beta function in other schemes. 
To calculate the beta functions we perform 



9_ 



d /87r2\ dp A 
dp\gi J dlog{fi/A) 

d /87r2\ dp ^ / /A\ A\ , , 



^These are strictly correct in the N = 2 examples and the KW fixed point. We adopt 
the definition here to get a handle on the non-SUSY dynamics. 
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We have reinstated the factor of Nc in the expansions. With a naive use 
of the NSVZ expression for the Wilsonian beta functions one may have 
interpreted this result for /3_ as the SUSY breaking parameter W20 changing 
slightly the value of the anomalous dimensions 7a,_b ~ — ^ + O ^^^20^^ . 
But this is not matching with the analogous calculation for /3_|_. Hence this 
solution does not respect the NSVZ expression (as expected). Also, notice 
that while in the SUSY case, the beta functions receive corrections O ) 
we have here an example where the SUSY breaking parameters produce 
lower order corrections O Let us move now to IR observables. 

6.1.5 K- Strings 

We will follow the treatment in [TO]. We need to evaluate the action for a 
D3 brane that extends on the manifold T, = [t,xi,9 = 9,f = 2Tr — if]. The 
D3 brane is sitting at /3 = but can move on the angle ^, so that it will 
minimize its energy. The (string frame) metric seen by such a D3 brane is 

ds^^a = ^ \dxl^ + nJiJ^ [dx^ + sin^ x {dO^ + sin^ Od^^)] \ (6.22) 



ho 

where we have written ip = 2x + t^, and used the values of the functions at 
p = 0: 

= = e^m = Q^ $(o) = ,/>o, a(o) = i. (6.23) 

We have additionally written ho = h{0) = 1 - e^'^^o-Z'J'oo^ 
The RR field and its potential are, 

Fsls = 27V,sin2 x^2^dx, Csjs = (x - ^-^] ^2, 



n2 = sine d9 A dip. (6.24) 



Using eqs.( |6lH6^ and the fact that 6'(0) = $'(0) = we find that the NS 
potential B2 vanishes. 

We will turn on an electric field F2 = Ftx dt A dx in the space-time 
directions. Then the Born-Infeld-Wess-Zumino action gives an effective one 
dimensional lagrangian, 



(6.25) 
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This is equivalent to eq. (9.8) of ^9j, with modifications which result from 
the U-duality, 



hi 

2 



hi 
2 



,24.(0) 



g20o 
hn 



(6.26) 



The rest of the discussion then proceeds as in pj] . We impose the equation 
of motion for Ftx and quantize it to be an integer multiple of the tension 



of the fundamental string, = -2 • The resulting tension follows an 

approximate sine-law, as in the whole baryonic branch, including the KS 
solution. This also happens for D5 solutions in section 8 of reference |12j . 
The influence from the SUSY breaking parameters enters only through 



kTf 



the modifications (6.26). 



6.1.6 The Non-SUSY Seiberg-like duality. 

We will follow the treatment in the SUSY case, as developed in pO|. The 
basic idea is go back to the quantity bo{ip), computed as specified around 



eq.(6.6) and compare with what occurs in the SUSY case. The Seiberg 
duality is identified with a large gauge transformation such that 60 — )■ 60 =t 1 
and the charge of D3 branes changes by ±Nc- 



Consider the Page charge od Section 6.1.3, a large gauge transformation 
on B2 will change 60 in one unit. This translates in the change of Nc units 
in the Page charge. This works exactly as in [20]. 

Let us now study how the Maxwell charge 'sees' the Seiberg duality. We 
will focus on the UV part of the background, where the cascade is known to 
work in the SUSY case. Following the steps described in Appendix [Gj we 
have 



bo 



TT 



626^'' - 64(0 + cosV'o)e'' 



Nr 



IT 



[{f + k) + {k 



with (using the explicit values for 62,64 
e*/2/ji/2 



/ 



2iVc 
e*/2/iV2 



60 



2N, 



6'(6 + cos V'o) 



o2$ 



-[6'(6-l) 



,2$ 



8 



iV2 

4 



TlNr 



' /) COS V-o] 
(6.27) 

g2g+2h^/j^ 
g2g+2/i^/j^ 

(6.28) 
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Now, it is interesting to notice that — far in the UV — the Maxwell charge 

QMa.,DZ = ^e^^'+^/^+^^cD' = ^^^^h\h + cos Vo) - N,h, (6.29) 
vr 47r 

changes under a change in 60 as, 

60 ~ 60 ± 1 QMax,D3 ~ QMax,D3 =F ^c- (6.30) 

Specially, notice that for large values of p the 'correction-term' 6'(6 + cos "00) 
is quite suppressed. This 'correction' is more suppressed in the SUSY case, 
where b' ~ e~'^P, in contrast to our non-SUSY solutions, where b' ~ e~^^/^. 
The 'Seiberg duality', associated with a large gauge transformation of index 
k that changes the Maxwell charge in kNc units is better approximated in 
the SUSY than in the non-SUSY case. Nevertheless, in both cases, the 
transformation is good at leading order. 

So, as expected, far in the UV we could think that the decrease in the 
Maxwell charge is interpreted as a non-SUSY version of Seiberg duality that 
is at work here. 



6.1.7 Domain Walls 



Let us compute the tension of a domain wall as the effective tension of a 
five brane that sits at p = and is extended along Sg = [t,xi,X2,9,'f,'4']- 
Before the U-duality for field theories of type II, we use the background in 
eq.( |2.6[ ) and obtain that the induced metric on such five brane is (in string 
frame) 



„<S> 



=23 



,2fe 



"^) + — (d'0 + cos0'd(^)^ 
The induced tension on the three dimensional wall is 



27r2rz55e2*+2s+^' 



p=0 



V2 ■ 



(6.31) 



(6.32) 



which is unchanged from the SUSY result. 

After the U-duality, in the background of eq.(2.12), we place a similar 
five brane, the induced metric is, 

1 



ds 



ind 



dxi 2 + y h 



2g 2k 



— {de'^ + sin^ ed^'^) + ^{di; + cos Odipf 



(6.33) 
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There is also an induced B2 field, 

B2 = ^Vhe^^+'^l%z{p) sin 9de A d(f. (6.34) 

In order to have a gauge invariant Born-Infeld Action, we must add the F2 
field on the world-volume of the brane. Indeed, the change due to a gauge 
transformation of the B2 field is cancelled by a (non-gauge)-transformation 
onF^l 

B2^B2 + dAi, F2^F2- dAi. (6.35) 

Hence, we need to turn on the worldvolume of the brane a gauge field 
strength, 

F^- = -'^e'^+''/%ip)sme. (6.36) 
This implies that the BIWZ action will be 

5 = -rTO(47r)2^— d^+^x (6.37) 



which gives the same effective tension as in eq.(6.32) and the same as in the 
SUSY case[^ Then the tension before and after the U-duality is the same. 
As a side remark, one may wonder if it is possible to fix the value of 63 at 
/9 = using some physical criterium. Though it is not an invariant quantity, 
the small p expansion of 

B^.uB^'' ~ + • • • (6.39) 

suggests that we should take 63(0) = as in the SUSY case. 



^One can also add a field strength F2 such that aside from cancelling the gauge- variance 
of B2 adds a kind of 'magnetic charge' to the domain wall or a Maxwell-like term in the 
Minkowski directions. We will not consider the addition of these extra components of F2 
as they will typically raise the energy of the wall. 



Notice, that in the SUSY case we have (using eq. F.l I 

h-j, = ~Ke^'^^'^h~^^^ cosa, (6.38) 

which vanishes for p = 0. 
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6.1.8 Central charge 

We will calculate the central charge of this non-susy solution. We should 
follow the usual procedure of [21], that requires a reduction to five dimen- 
sions. However, an equivalent treatment presented in [22] indicates that for 
any string-frame metric of the form 

ds^ ~ aip)dxla + a{p)l3{p)dp^ + gij (p, y)dy'dy^ (6.40) 



we define 



Vint = I dy^/deUHj, H = e-^^a^.^^- (6-41) 



and the central charge (for d = 3) is given by c ~ ^ [h^)'-^ ' ™ ^^'^ ^^^^ 



{2h' + 2g' + 2$' + A:' + 2^ 



(6.42) 



In the IR, the explicit expansion for the central charge is 

1 ....... /^2<i.^ _ ^2^,-^ ^2 L 2i>^ (_i6 _ l5/,^fc2 + 12/if) 



9 

+ e^'t"'' (28 + I5hik2 - I2hj + 9v 
and in the UV we have 



p' + 0{p'), (6.43) 



It is immediately clear that the SUSY-breaking parameters have no effect 
at the leading order in the UV. However, in the IR the question is more 
subtle. Although none of the explicit SUSY-breaking parameters appear in 
the leading term there is an effect. This is because in the SUSY case there 
are only two independent parameters, so that fixing hi and (po is sufficient 
to determine <I>oo- In the non-SUSY case the discussion of sections [3 . 1 1 and [l] 



suggests that there is one more parameter, which breaks SUSY. This means 
that even with fixed hi and (f)o we can expect that <l>oo varies as a function 
of the SUSY-breaking parameter. Indeed, when in appendix |D] we compare 
SUSY and non-SUSY numerical solutions with the same hi and (po, we find 
that $oo changes. 
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6.1.9 Force on a probe DS-brane 



We will now consider a D3 probe brane that extends in the Minkowski 
directions and is free to move in the radial direction as suggested in [16] , 

: [t,xi,X2,X3], p{t). (6.45) 

the induced metric and RR four form field are obtained from the string 
frame version of eq.(2.12), 

„2 



^2* 



-K^^dt A dxi A dx2 A dx^. 
h 



this gives an action for the D3 bran^ 
Sbiwz = -Tb3^3 / dt 



o2* 



1 — he^^n''^ — K^r- 

h ' h 



(6.46) 



(6.47) 



We then approximate this for small velocities and change to the variable 

dr = e^'^^^'^dp and get 



(6.48) 



S = TdsVs I dt 
the force on this probe is then 



1 + Ke 



^ (1 + Ke*)2 ■ 
In the IR, the explicit expansion for this force is 



(6.49) 



/ 



2y/2 



'to. 
e 2 



+2$o 



(4 + 2,vl) 



3 (e<^o + e*-)^/ii/^ 



P + 0{p^), 



(6.50) 



and in the UV we have 



3 e 2 



(9e4*- + 4d 



301 



2 5/2 



8\/6c 



5/2 



e-i0p/3 + o(g-i4p/3)^ (6.51) 



^Notice that in the way things have been defined, the action for a D3 has the WZ term 
with the same sign as the BI term. See eq.(2.13) in the paper [9]. 
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As expected, the force vanishes quickly in the far UV, where the solution 
approaches the KS background. Also, notice that in the radial coordinate 
r ~ e"^''/^, the force is / ~ as obtained in [23j. The SUSY breaking 
parameters do not influence this small force. In the other hand, the breaking 
of SUSY explicitly changes the value of the force in the IR, as expected. 

6.2 Field theory comments 

This section relies on the ideas of [24j-|25j. but most fundamentally on the 
analysis of the paper [25]. Similar ideas that may be useful in thinking 
about our string backgrounds have been put forward for example in j26j . 
This paper studies non-SUSY deformations of A/" = 1 SQCD. We use this 
to analyze the quiver field theory of type I. This is as we discussed, a non- 
SUSY deformation of the KS-quiver. In the SUSY case, the KS-field theory 
can be understood as N=l SQCD with gauged flavor group and a quartic 
superpotential (see for example [28]) and due to this, the results of [25] are 
important to us. The qualitative results of the paper [2H] become quantita- 
tively accurate once we take the SUSY breaking parameters much smaller 
than the relevant scale of the problem, namely AgQgJ^ In our case, this 
is reflected in the smallness of the coefficient W20 ■ 

In this case, lots of the structure of Seiberg's SQCD [29] remains. Par- 
ticularly interesting to us is the fact that for SU (Nc) SQCD with Nf flavors 
and with Nf = Nc, there exists a vacuum which breaks spontaneously the 
C/(l)-baryonic symmetry and this vaccum persists in the non-SUSY analysis 
of |25j . This will be relevant for us as the case Nj = Nc is associated in the 
SUSY case with the last step of the cascade. We then argue that our ge- 
ometry describes a situation where SUSY is broken by gaugino masses and 
other VEV's and the baryonic symmetry is broken by the vacuum state. 

In a bit more detail, the authors of [25], added to the SQCD lagrangian 
a term of the form 



^ Ofer Aharony explained us that the soft breaking mass terms for squarks could have 
different signs under a Seiberg Duality, see [27] • This technical subtlety seems to play no 
role in our analysis 



L = LsQCD + AL, 



AL 




(feMQ{Q^e^Q + Q^e-^Q)+ (fOMgS, (6.52) 
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where S is the superfield S = Tr[WaW°'], Mq is a vector multiplet whose 
D-component equals the mass of the squarks (— m^) and Mg is a chiral 
multiplet whose F-component is the mass of the gluino. The authors of [25j 
argued that to leading order in the SUSY breaking parameters Mq , Mg one 
can write an effective lagrangian in terms of mesons M, baryons (B, B) and 

AL^ J (feBMMQtr[M^M]+BbMQ{B^B+&B)+ j (feMgS+.... (6.53) 

The idea is then that one should supplement the usual actions and superpo- 
tentials discussed in the SUSY case with the SUSY breaking terms above. 
In particular, in the case Nf = Nc we will need to minimize the potential 
term coming from eq.(6.53) together with the potential coming from the 
SUSY superpotential 

W = Wtree + Wguant = KTr{M^ M) + ^(det M - BB - A2^^). (6.54) 

Therefore, the vacua of the theory are those that minimize the potential 
coming from the tree level superpotential, together with that from the SUSY 
breaking term, all subject to the constraint in Wguant- The result is that 
in the non-SUSY case, one finds one vacuum state where the baryons get a 
VEV and the mesons are at the origin of the moduli space, M = 0. 

In this way, we have argued that our solution, which breaks SUSY due 
to masses for the gauginos has very similar behavior to the KS-cascade (ac- 
tually to the baryonic branch in [3j). We found that many non-perturbative 
aspects behave very similarly as the SUSY case: the expression of the do- 
main wall tension is basically the same as in the SUSY case. Of course, 
numbers will differ as the functions in the IR pick the influence of the SUSY 
breaking terms. The tensions for k-strings gives an approximate sine-law, 
again with the SUSY breaking entering the value of the tension. In the UV, 
the beta function for the gauge couplings of the quiver and the leading order 
of the central charge behave at leading order in the UV like their SUSY coun- 
terpart, but in the case of the beta functions, the first correction is purely 
coming from SUSY breaking contributions. The Seiberg duality (identified 
here with the change of the Maxwell charge under large gauge transforma- 
tions of the NS B-field) behaves very approximately as in the SUSY case. 
One can probably make an argument for self-similarity as presented in |28j . 
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Obviously, what happens is that the SUSY breaking terms, for example 
the gaugino mass indicated by the quantity W2o, are not important at high 
energies. They enter some IR observables, correcting but not changing the 
qualitative behavior expected from the SUSY example. This suggests that 
we need to think that our SUSY breaking scales are smaller than our strong 
coupling scale. Hence, the phenomena are the same as in the SUSY models, 
but numerically there will be differences. All this is in line with the analysis 
of [2l]-[25]. 

7 Conclusions 

In this paper we have used analytical (UV and IR series expansions) and 
numerical methods to construct smooth backgrounds dual to particular non- 
SUSY field theories. The field theories in question can be thought of as 
softly-broken-SUSY versions of the field theories appearing in twisted D5 
branes and Klebanov-Strassler quivers. 

We presented some details of the derivations, involving a U-duality, a 
careful numerical procedure and a detailed study of many observables at 
low and high energies. All this supports the field theory interpretation 
discussed towards the end of Section [6j In other words, the dynamics is 
basically the SUSY one, but with interesting details and deviations coming 
from the soft-breaking terms. 

Various things come to mind that would be nice to study using the 
backgrounds presented here. Before the U-duality, it would be nice to study 
the effect on k-string tensions, domain walls and the confining behavior 
of the Wilson loop, as there exist in the bibliography various results for 
Af = I Super- Yang-Mills being softly broken. Also, it would be nice to 
study the effect of the mass terms responsible for the breaking on the glueball 
spectrum. It would also be interesting to see if our solutions can be of any 
help for the interesting line of metastable broken SUSY, in the sense of 
providing a good set of UV boundary conditions that break SUSY. This 
may be used to get ideas on the singular IR behaviors obtained in [33] , [23] . 

It would be interesting to calculate the mass spectrum and compare it 
with the result of the analogous glueballs in the KS/baryonic branch solu- 
tion [30]. Also, it would be nice to find numerically the expression for the 
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massless glueball corresponding to the spontaneous breaking of the bary- 
onic symmetry. The spectrum of mesons is also of interest. In particular, 
comparing the masses of the lightest scalar meson and the lightest vector 
meson we could learn if the non-SUSY background presented here provides 
a plausible holographic dual of nuclear forces [31] , [32] . 

Another problem that we are not addressing here: it is known that 
taking the integration constant c+ — )■ oo leads, in the SUSY case, to the 
Klebanov-Strassler background (see [9] ) • It is interesting to see this working 
numerically and to compare the solutions found here — in the limit — with 
those found in the past by first order fluctuations of the KS system (see 
|34j ) . In this line, a nice problem would be to find the recent solution by 
Dymarsky and Kuperstein [35] as a scaling of ours. 
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A Appendix: Technical aspects of the SUSY back- 
ground 

We write in this appendix various technical aspects of the supersymmetric 
backgrounds. As explained in Section [2] one changes the basis of functions 
from ^,h,g,k,a,b into P,Q,Y,^,T,a in order to decouple the non-linear 
system of BPS equations. As explained in |13| . the change of basis functions 
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IS 



4e'^= ^\ e^f = Pcoshr-Q, e^'' = AY, 

P cosh T — Q 

P sinh r , ^ , 



(A.l) 



P cosh T — Q 

Using the relations above, one can solve for the decoupled BPS equations, 



Q(p) = (Qo + Nc) coshr + iVc(2/0COshT - 1), 



sinh t{p) 



1 



sinh(2/) — 2po) 



, coshT(/9) = coth(2/9 — 2pq), 



4$ ^ e"'^° cosh(2po)2 



(7 = tanhT(g + A^c 



{2N,p + Qo + N^ 
sinh(2p — 2po) 



(A.2) 



and the master equation (2.8). Solving the master equation in the UV (2.9) 
and plugging back into eqs.(A.l)-(A.2) the background functions read at 
large p, 



^2h 

e ~ 



=2g 



c+e 



4p/3 M 



A^2e-4p/3 

4 +TP''-^) + ^ 



(16p2-l6p+13) + 



-8p/3 g 

— (c_-c+(2 + |)) 



4 

„2fc 



f(2.-l) + -^(16p^-16p + 13) + — 

^ - ^^^(4p - 5)^ + ^(c,(^ - 1) - c_) 
6 24c+ ^ ^ ^ 3 ^ +^ 3 ^ ^ 



-(c- + c+(2-^))_ 







3^2 -8p/3 



+ 



3^4g-16p/3 
5124 



(2048/)^ + 1152p2 + 2352p - 775) 



2p 



b = — ~ Ape-P + Ape-^P 
sinh(2p) 



(A.3) 
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The geometry in eq. (2.6) asymptotes to the conifold after using the expan- 



sions above. In the IR we have using eq. (2.10) and (A.l) 



2h 



=23 



45 



-6/., + 15iV. ^^^^^^ 



Inj^ ^ 1575/if 






hi {hi - 
8 


4 


g4(<&-4.o) _ 


64iV> 


a ~ 


1+ (-2 + 



10/ii 



+ 



315/if 



405/if 



(A.4) 



8iV, 
3/n 



2 , 2 (75/1? - 232/if iVe + imhiNj+^m^ ^ ^^^^^ 



45/if 



2p 



smh(2p) 3 45 



(A.5) 



This space is free of singularities as can be checked by computing invariants. 

B Appendix: Euler-Lagrange equations of motion 

Here we write the full equations of motion, for reference. We start with 
the effective Lagrangian and the constraint and then write the equations of 
motion. We set A'^c = 1 for simphcity. 

The effective Lagrangian is L = T — U, with 

T = --l-e^^je^^ {af + {bf - 8e2(^+^) hg' {2h' + k' + 2$') + {gf 
128 L L 

+ 2h' {k' + 2¥) + {h'f + 2$' {k' + }, (B.l) 

256 L V ; " \ " 

+ e'^9]y2 _^ 4g2/»^2 _ gg2(g+/i+fc) _^ ^^4g+2h _ ^2g+4k _^ ^^2h+ik^ 
- Aabe^ml { + ^e^h^ + ^^2^2^2{g+h) ^ ^4g^2 ^ iQ^4h^2 

(B.2) 



- IQe^C^a+h+k) _ Q^^2{g+2h+k) _|_ g4(g+fc) _^ IQ^^h+k) 
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The constraint is 



= T + U 



-2(s+/i-*) 



The equations of motion are: 



(B.3) 
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1 

-e 



-Ag-2h 



16e^»+2'» {g'f + ^2e^9+2h+2k _ iQ^2h+Ak _ ^Q^2h 



(B.4) 



h" 



-2g-ih 



(a')2 ^^9+2h + ^4g2s+4fc ^ ^4^23 _ 4„3j,g2p ^ ^^2^2^25 



■2h (iJ\2 



- 8a^e^9+2h+2k ^ 4„2g4g+2/. _ 2a^e^<^+^^ + 2a''e''3 

+ 4a2e2''+^fc + 4a2e2'' - 4a6e2f - 8a6e2'* + 

_ gg2s+2/x+2fc _^ g2g+4fe _|_ g2g^ 



(B.5) 



- 8a2e25+2'^ _ 8G2e65+2'^ _ 2a2e^^+^'= - 20^6^^ + lQabe^3+''^ 



- l&e'^9+i^k'<t' + e^3+^^ - + lQe^^+^^ - 16e 



Ah 



(B.6) 



-Ag-Ah 



a^e^9 _ 4„3^g4, ^ 4^2^,2^49 ^ 8a2e2^+2'^ - 16a6e2^+2'^ 



+ 2a2e^5 - 4a6e45 + 2 {b')\''s+'"' + 8b''e^s+^'' - X^e^^^^^g'^' 



(B.7) 
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+ 4ae2'*+^'= + 4ae2^ - be^^ - Ahe^^^ (B.8) 

C Appendix: Explicit expansion of the functions 



Here we include the explicit solutions for the expansions (3.1) and (3.7). In 
this section we again set Nc = 'i. 

C.l uv 



p29 



„2h 



2k 



4$ 



c+e 



4p/3 



- ^{-12^^11 [(32P - 6)Qo - c+T^|o(8p + 93)] - 12c+VF|oQo 
+ 724^206^''° + 12Qc\W^qP - 2Qc\W^Q + 124$3oe-^*°° - 72p 

- 2AHl^ (32p2 _ I2p + 15) - A^QI + ^^e'^'''^ + 0(e"^^/^) (C.l) 

- [c+VF|o(88/J + 75) + (32p - 6)Qo] - 396c+PF|oQo 

+ 26441^206^''° + UOclwioP - 626clWio + 124«>3oe-^*°° - 72p 

- 2AHf^ (32p2 _ Up + 15) - 48Q2 + gje-^^/^ + ©(e-^^/^^ (C.2) 

+ ^ - ^ [4//n(16p - 9) (3c+< + 2Qo) + 16Q^ 
+ 84c+VF|oQo - 724^206^^° - UOcIW^qP + 1904^20 - 24p - 9 
+ 44«>3oe-^*- + 8Hf, (32p2 _ 36/? + 27)1 e~*P/'^ + 0(6"^^/=^) 

(C.3) 



e'*- + ( <I>30 



-8p/3 



, ^lo (324^^30 - 3(64p + 25)e^^-) ^/ -l6p/3^ 

+ 24c2 e +cy(^e ; 



(C.4) 
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1 



48c2 



4c+<^ i^u [96pe2^° + Wio (iGOp^ + 552p + 495)] + 2Qo \l2e^''° 



+ Wio{20p + 87)1 [ + 4 [^2o(391 - 480p) - 288W^|oe^''°] 



+ 12W^20 



2ffii(40p + 27)Qo 



+ 6Hj^ (32p2 + 36p + 43) + 8Ql - isj e-^^^/'^ + 0(6"^'') 



(C.5) 



6 = ^T^soe-'"/^ + 



4e^- + M)(20p-23)-^^^ 
6 6c+ 



-2p 



I 3W2oe- - J 
^ 512c2 ^ 



-3456iJii (3c+W|o + 2Qo) - 3240c+W|o<3o 



+ 19241/40^^20 + 2819clwio - 35712i?i2i - 576g^ + 882 



16pe^*°° [l44ifii {3c+Wio + 2Qo) + 24c+W|o<3o 

4^20 (486^''° + 269^2^0) + 1728iJfi - 18] + 30244^^206^''°+^*°° 

128p2e4*- (721^2^ - 54W^2'o) - 484$3o|e-'°''/^ + ©(e"^^) (C.6) 



We can see the effect of the SUSY-breaking parameters by looking at 
functions of the form A{e^^) = e^^ — e^^^usY^ where g corresponds to the 
full solution and ^suSY corresponds to the SUSY case with Qg = —Nc and 

Po = 0. Note that in general only one of the two solutions will go to the 
regular IR if we start with a SUSY solution and turn on one of the 
SUSY-breaking parameters while keeping e.g. c_|- fixed, we will have to 
change c_ to recover the regular IR. Those SUSY-breaking parameters which 
have non-zero values in the SUSY case are expressed here in terms of e.g. 
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A(e2^) 







^ 24c+ 








Ai7n((32p 



-2AQo + 2Aii"ii(8p + 93) + 8p + 95 

|o(60p-13)e^*°°+6A$3o] 
Ai^ii ((32p - 6)AQo + 64p2 - 56/9 + 36) + AQo (4AQo + 16p - 11) 



+ AHfi (64/9^ - 24p + 30) 



.e-4p/3 + o(g-8p/3) 



(C.7) 



+ AHuP + 



4$o 



3c+W^2oe 66AQ0 + 2Ai7ii(88/9 + 75) + 88p + 9 

6A$3ol 



V 4 96c, 
- 4 [1321^206^"°+^*°° + W^lo(220p - 313)e^* 
+ 6e^*- 1 Ai^ii [(32p - 6) AQo + 64^^ - 56p + 36] + AQo (AAQo + 16p - 11) 

+ AH^^ (64p2 - 24p + 30) 1^ e'^''/^ ^ o(e-8p/3^ 



1 



3^+^- + T2c7 



(C.8) 

14AQo + 2Aii'ii(16p - 9) + 16p - 23 



+ 4 [36W2oe'''<'+^*- + 5Wio{l2p - 19)e^*- - 2A$3o] 
- 2e^*- {2Ai?ii [(16p - 9) AQo + 4 (Sp^ - I3p + 9)] 
+ AQo (4AQo + 16p - 17) 

+ AHf^ (64p2 - 72p + 54) }^ e"^''/^ + 0{e-^P^^) (C.9) 

T/t/2 

A(e^*) = A$3oe-'''/' + ^ [324A$30 - 3(64p + 17)6^*-] e-^^ + Oie-'^"/^) 

(C.IO) 
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Aa = VFsoe"^^/^ + — 



-c+ (3e2^° + 5W|oP - 3) + 2W20 {2AHu + 1) p 



-2p 



+ 



48c2 



4(20p + 87)AQo + 2Affii (l60p^ + 552p + 495) 
+ 160p2 + 472p + 147 + 48 [AAHupe^P" + AQoe'^^" + (2p - 1) (e^''" - l)] 

+ 4 [^lo(391 - 480p) - 288Ty|oe^^°] + 6H^2o|+(80p + 22) AQ^ 
+ AAHn [(40p + 27)AQo + 96p^ + 68p + 102] + IGAQl + 12AHf^ (32p2 + 36p + 43) 
+ 96p2 + 28p + 6l|^ e-^O''/^ + o{e-^'^P/^) (C.ll) 

Ab = \w-2ae-'^Pl^ + [c+ {p [24 (e^^" - l) + W|o(20p - 23)] + 6^40} 



-2p 



- 12T^20P (AQo - 1)J + ^^'°^^2 (^4 |48W2oe4*- [(16p + 63)e2''<' + 4^40] 

+ W2^o (640p^ + 4304p + 2819) e^*°° - 48A$3o| 

- 24c+M/2V^*°° [(16p + 135) AQo + 144Ai?ii(2p + 3) + 128p + 81] 

- 186*^*^- |64Ai?ii [(4p + 6)AQ<, + 8p2 + 20p + 25] 

+ 128(p + 1)AQ<, + 32AQ2 + 64Ai72^ (Sp^ + 24p + 31) 

+ 128p2 + 240p + 289| j e-^^^l^ + 0{e-^^Pl'^) (C.12) 



C.2 IR 



,25 ^ ^ _ - 



^ - g [4A;2 - /ii + 4) + ^ (^i + 4) J p2 + 2QjgQ^ 1 1600/^1^2 (3^1 + 8) 
- 8/i? [210A;i + 144i;2 (w2 - 3) u;2 + 3v| (l05u;| + 168u;2 + 580) 

+ 4 (55ii;i + 360^2 + 652) - mh\k2 {wl + 18t(;2 + 18) 

+ 3/if (75«;^ + 432w;f + 488u;i + 960it;2 + 1520) 

+ 16 (405v| + 1592t;f + 656) Ip^ + O (p^) (C.13) 
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e = — p -I 

2 ^ 72 

1 



24^2 - 3/ii {?>wl + 4) - T-(9ui - 4) 

"•1 



24/it 



5250A;i + 432i;2u;2 (23u;2 + 57) 



„2k 



907200ft,f 

- 9i;| (I89u;i + 2016u;2 + 940) - 20940w| + 8640^2 + 5680 
+ 480^1 A;2 {Qvl - 172) - 360/i?A;2 (627u;i - 432u;2 + 44) 
+ 9/if (2007u;| - 6624«;| + 11352u;i - 5760u;2 + 1520) 

+ 16 (I9359i;| + 27288v| - 22672) jp^ + O [p^) (C.14) 

h\ , n 2 , / o „x 15 



y + k2P^ + ^Y^|5^i^2 (27^;^ + 2) + —hfki {9wl + 36w2 + 22) 
hj \l75k^ + Y1V2W2 {w2 + 4) - 30u| + 30w;^ + 120w;2 + 208 



16 



h\ {w\ + %w\ + 36u;| + 80iU2 + 80) 



+ 9v\ + Z^vl - 528|p^ + O (/) 



(C.15) 



4e4</'o 



,4* ^ ^4,0 + (3^2 + 4) p^ + ^ I -60/.ifc2 (3^1 - 8) 

+ 3/if 3i;| (9w^ + 36w2 + 40) - 36?;2w;2 («^2 + 4) - 176 
+ 4 (243i;| + 672uf + 944) + O (/) 



a = 1 + 'u;2p + 



W2 



150/11 A;2 - 3/i? (6u;i - 9u;2 + 28) + 400 



+ 36v^ {W2 + 2) - 72^;2 {w2 + 2)|p^ + O (p^) 



(C.16) 



(C.17) 



6=1 + i;2p^ 



90/if 



t'2 



30hik2 - hi (9w;| + 36k;2 + 60) + 176 



+ Qhlw2 {W2 + 4) + 72vl \p^ + (p^) 



(C.18) 



The effect of the SUSY-breaking parameters can be seen from the dif- 
ferences: 
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1 



24/ii L 
+ 



-4/ii (3AA;2 - AAw2) + Shj {Aw2 - 4) Aw2 + 4 (4 - SAus) Av2 



, . -24/i? [320A/C2 {Aw2 + 7) + 210AA;i - 4Av2 (69Au;| + 224) 
60480ft.f [ L 

+ 3Avl (l05Au;| - 252Aw2 + 584) + 8 (112 - 129Aw2) Aw2 

+ 16A?;2 (l215Av| - 3240Ai;| + 3216Ai;2 - 2944) 

+ 9hiAw2 {75Aw^ - 168Aw^ - 368Au;2 + 896) 

- 96/t? [3AA;2 (SAu;^ + 70Au;2 - 56) + Au;2 {-75Awl + 126Au;2 + 184) 



+ 64/11 



3Afc2 {75Av^ - 100At;2 + 264) - 126Aw| {5Aw2 - 6) 



+ 552Ai;2AiU2 + 464Au;2 

1 



(C.19) 



24/ii L 
+ 



8hi {Ak2 - 2Aw2) - Shj {Aw2 - 4) Aw2 + 4 (4 - 3Au2) Av2 
1 

302400/if 

+ At>| {-567 Awl - 3780Aii;2 + 7032) + 4At;2 (l017Au;| - 3136) 



2Ahl -80AA;2 (209Au;2 - 490) + 1750AA;i 



+ 40Au;2 (157Au;2 - 1288) 



+ 32/11 



5AA;2 {9Avl - 12At;2 - 4352) 



- 4 (l89Av| (3Au;2 + 10) - 4068Au2Au;2 + 856Au;2) 

- 24/i? [5AA;2 (627Au;i - 2940Aw2 + 3136) 

- 4Au'2 (669Awi - 5670Ait;2 + 14872) 

+ 9hlAw2 (669Au;| - 7560Aw| + 29744Aii;2 - 49280) 
+ 48Av2 (2151Av| - 5736At;| + 6704Ai;2 + 7936) + O (p^) 

(C.20) 
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A(e 



2k\ 



Ak2p^ + 



7560/if 



15AA;2 {'iAwl + 14) - SAw2 {Awl - 6) 



+ 12hl (l20A/c2Att;2 + ITSAfc^ + 12Au2Au'i + QAvl - bQAv2 - 30Aw|) 
+ IQhi 15AA;2 (27At;| - 36At;2 - 106) + 16 (18Ai;2 - 29) A'u;2 
+ 16At;2 (27A?;| - 72Aw| - 468Ai;2 + 1072) 

- 27h\Awl {Awl - 12) |/ + O {p^) (C.21) 

400 



4e^'^''At>2 , , ,,2 e 
(3At;2 - 4) p2 + 



2,hl 



135/i| 



4/n 



AA;2 (-45Af| + 60Au2 + 100) 



+ 36 (3A^;i - 8At;2 + 4) Au'2 + 3/i? 3At;^ {9 Awl - 4) 
- 8Au2 {^Awl - 20) + 36Aii;| + 4Ai;2(243At;| - 648At;| 

+ 1536AV2 - 1664) + O (p^) (C.22) 



Aa = Aw2p^ + 



90/if 



4/n 



100AA;2 + {^Avl - 30At;2 - 40) Au'2 



+ 6/i? 25AA;2 (Au;2 - 2) - 24 (Au'2 - 5) Au;2 + 32A^;2 (3Au2 - 10) 



Ah = Av2p + 



3^?Au2 (6Au;i - 45Au;2 + 116) }>/ + O {p^) 
1 



(C.23) 



90/if 



Afc2 (20 - 30Ai;2) + 16 (3Au2 - 5) Au2 
+ 2>hl Av2 {3Awl + 4) - SAu; 
- 8Au2 (9At;| - 18Ai;2 + 20) Ip'' + O (/) 



(C.24) 



D Appendix: Details of the numerical analysis 

Here we shall discuss in more detail our approach to connecting the given 
IR and UV expansions numerically. We start by noting that we have chosen 



to solve the equations of motion (B.4-B.9) starting from the IR boundary 



conditions. As the IR parameter space is much smaller than that of the 
UV, this makes a search for solutions with the correct behaviour much less 
computationally expensive than if we started from the UV. 



We use as our boundary conditions the IR expansions (C.13-C.24), ex- 
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tended up to order . Using NDSolve in Mathematica 7 we then are able 
to generate numerical solutions which extend into the UV. We start at 
piR = 10~^ as we found that in the SUSY case this gives approximately 
optimal accuracy. We use 40-digit WorkingPrecision in NDSolve. 

Comparing the numerical solutions obtained by this method, with the 
known behaviour in the SUSY case, suggests that the results are trustable 
up to /) ~ 11. This is supported by the observation that the constraint 



(B.3) is almost completely satisfied over this range. More explicitly, we find 
T + [/ < 10^* throughout this range. In fact it appears that the numerical 
solutions only fail when h decreases past ~ 10^^. In the SUSY case (in 
which h ~ e"^'') this does correspond to p ~ 11, but in the non-SUSY case 
(with b ~ e^^^^) it occurs further into the UV. 

In the IR we have five parameters {hi, (po, W2, k2, V2} which we can ma- 
nipulate although we set (po = (along with Nc = 1) without loss of gener- 
ality. Given a value of hi we want to study the effects of the SUSY-breaking 
deformations for the remaining three {^2, ^2, ^2}- We find that for a general 
deformation of these IR parameters the resulting solution does not exhibit 
the expected UV behaviour. Initially we find that the general behaviour of 
solutions in this parameter space has 

br^±e'^P and e^^^ ~ e^^ ~ e^^^ ~ e^''/^ p^^) 

going into the UV. The e^''/^ behaviour appears to be suppressed by a very 
small numerical factor relative to the expected e'*^/^ term, and in fact is 
not visible in plots of g, h and k themselves. It is apparent, however, if we 
examine quantities of the form 

g2fe _g2fcsuSY _ g8p/3^ (D.2) 

in which the e^''/^ behaviour (almost) cancels. 

Given a value for one of the three non-SUSY deformations, we believe it 



is possible to obtain the desired UV behaviour (C.1-C.6), i.e. 



with the correct choice of the remaining two. In practice it seems easier to 
vary three but keep one very close to its starting value. 
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Having obtained a numerical solution with the correct UV behaviour, 
we look to determine the corresponding values of the expansion coefficients 
in the UV, i.e. 

{c+, c- , ^>oo, Qo, Po, i^ii, VF20, ^-so, Vm}- (D.4) 
We define the mismatch function 

E [//"""''(Pmatch) - /r^^°^'°"(tach)] ' , (D.5) 



m 



with fi G {g, h, k, <I>, a, b, g' , h' , k' , <!>', a' , b'}. We then minimise m to match 
our numerical solution and a UV expansion using NMinimize (with 60-digit 
WorkingPrecision) at a large p value, Pmatch- 

With this setup and given the SUSY IR, NMinimize recovers the SUSY 
values for the UV parameters with an acceptable accuracy, even allowing all 
nine parameters to vary. The only restrictions we apply to the parameter 
space are c+ > and $00 > (/"o = 0. 

We now present a non-SUSY solution found using the above methods 
for one set of values of the IR parameters. It has the expected behaviour for 
all functions at least up to p ~ 11 (where the corresponding SUSY solution 
fails) and possibly as far as /) ~ 17. We first choose hi = 5 (and have set 
00 = as mentioned above). The corresponding SUSY solution has 

8 n 22 , 2 , . ^, 42 2 
"^=3/.r"' = "l5' ^^=5/^(^1-') = 25' "^ = -3- 

This results in an NMinimize output (with Pmatch = 6) of 

c+ « 1.6, c_ « 2.0 X 10^, 



$00 - 


a 0.076, 


W20 ^ 


a 6.9 X 10 


$30 ^ 


a 0.38. 



Qo ~ -1.0, po ~ -6.8 X 10-11, 

V40 w 2.7 X 10-^ Hii w 0.50, 

The associated mismatch value is m < 10-29. We take this as a good value 
for the mismatch as we know that the SUSY solution does indeed exist, and 



these values are in good agreement with (3.5). 

To obtain a non-SUSY deformation, we follow the procedure described 
and modify the three IR parameters {^2) ^^2) ""^2} away from their SUSY val- 
ues so as to manually scan the parameter space, until we gain a solution with 
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the correct UV behaviour. We find that a suitable choice of deformations 

ifO 



AA;2 ^ -2.471 x 10"^ Av2 ^ 2.574 x 10'^ Aw2 ~ 1.029 x 10"^. 

The minimization routine (again at Pmatch = 6) then finds that the UV 
parameters are modified from their SUSY values according to 

Ac+ -6.6 X 10"^, Ac_ pa 1.6, A<I>oo -4.0 x lO^^^, 

AQo ^ -1.5 X 10^"^, Apo ~ -7.1 X 10-^ AVFao ~ 5.2 x 10~^ 
Ay4o 5.6 X 10^^ AiJii 9.1 X 10-^ A^>3o « -5.0 x 10-^ 

again with a mismatch value of m < 10^^^. However, we are unsure of the 
precision of these values — they appear to be slightly sensitive to the value 
of Pmatch, and so should be interpreted with caution. We present plots of 
the functions (figure [T]) in the main text. 



E Appendix: Free Energy 

Consider the Euclidean action I for the wrapped D5 background of section 
([3|. The free energy is F = I/f3, where /3 is the period of the compactified 
time direction. 



2^ — -Z^grav ~l~ -Z^surf 
_ 1 

~ ~16^ 



- i ^KdJ^r- (E.l) 



M JM 



7W is a ten dimensional volume enclosed by a nine dimensional boundary 
Tir- The boundary T,r is taken to be surface at constant radial direction r. 
is the extrinsic curvature of the boundary, 

"^K = V^n^ = -^5^ iVsn^ = ^e-^/^e"'^ [9$' + 8(5' + h') + ik')] (E.2) 



The exact values used were Ak2 = -24 705 875 x 10"^^, 
A«2 = 25 744 091 286 331 971 640 358 x 10"^^ and Aw2 = 1 029 383 373 181 636 875 x 10" 
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where is the boundary outward normal vector, = ^J~(f^br ■ Using the 
equations of motion Xgrav reduces to a volume integral of a total derivative, 

= vohfi^ hm (ie2(*+5+^)cl>'') . (E.3) 

Explicitly, the surface term is 

Xsurf = -vohH^ hm I le2(*+s+'^) [9<1>' + 8(5' + U) + 4A:'] | . (E.4) 
87r r-!>oo 1^ 32 -'J 



Thus, 



2^ — 2grav ~l~ -^surf 



2567r 



lim |e2(*+9+'^) r8($' + 5' + /i') + 4A;'l j . (E.5) 



Equation (E.5) gives the value of the on-shell action in terms of the asymp- 
totic fields at infinity. It typically contains divergences and has to be reg- 
ularized. One way of doing this is to subtract the action of a reference 
background. In our case the natural choice is to subtract a supersymmetric 
background. We also require that both backgrounds induce the same metric 
at the boundary, S 



g 2 g^y^s = g 2 g^y^^ g 2 g^'^^is = g 2 ^^'''SU 

e 2 e =e 2 e e 2 =g 2 (E.6) 

and that the matter fields coincide at the boundary. In order to achieve the 
matching of the induced metrics and matter fields at the boundary we have 
to choose particular values for the integration constants of the supersym- 
metric background that we use as a regulator. We can then evaluate the 
free energy. 



vols 
'2567r 



_g2^.,+25.+2/^, ^8^/^ ^ ^ 8/i; + 4A:;) } . (E.7) 
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Using the UV expansion (3.1), to first order in W: 



'20, 



E 



vols 
2Att' 



20- 



(E.8) 



which agrees with the ADM calculation. A similar evaluation of the free en- 
ergy can be carried out for the backgrounds after the rotation. Due to the 
presence of and the Chern-Simons term the calculation is more involved 
and the equality of the energy before and after rotation cannot be expressed 



as simply as (5.12). Nevertheless, after plugging in the appropiate UV ex- 

cle'^P'>+^^^W2o as 



pansions we get, to first order in W20, Fi 



before 



after 



expected. 



F Appendix: Calculation of B-2 

In the SUSY case, we have 



Bo 



62 



/1I/2 L 



.P3 



cosa(e^^ + e''') - sina(e''" + e^^) 



-,12 



92 



with 



cos a 



cosh(2/9) — a 
sinh(2p) ' 



sma 



sinh(2p) 



(F.l) 



(F.2) 



This is not valid in the general non-SUSY case. We obtain the same 



as in the SUSY case (2.12), but the relationship to (F.l) requires the BPS 



equations, as does the consistency of the definitions (F.2). 



Instead, we must determine B2 by requiring that di?2 = ^^3- We assume 
it B2 has the same general structure as (F.l), 

B2 = hi{p)ef^ + b2{p)e'^ + h{p)e'^ + b,ip)e'^ + h{p)e^\ (F.3) 
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which results in 



(a63e^ + 2b5e^) ^^3 



(64-fe5)e-''-*cot^ 



/lV4 



/lV4 

-29-fe-f 



+ 



+ 



2/l5/4 

g ft 4 

2/l5/4 



,29 



536^^0' - 2a6ie-^?/ie2^ 

+ e'^ {A [64 (25' + 2/1' + + 26^] + 64/1'}') eP'^ 



+ 



2/l5/4 



g ft '« 4 
2/l5/4 



+ ^ . " ( &3e»/ia - 2a6ie-%e2^ 



+ {A [65 {2g' + 2h' + + 26'5] + 65/i'}) e^^^ 
-(64 + 65)e^/ia' - (a2 - l) fti/ie^'^-'' 



+ e' 



{h [62 (4/i' + + 2b'^] + e"''^. (F.4) 



Comparing with (2.12), we see that the e^^^ component of is zero, from 
which we immediately obtain that 64 = 65. The e''^^ and e'"^^ components 
of (F.4) are then identical, as are the e^^^ and e"^^^ components. This is also 
the case in ^^3, so we are left with four independent equations. 
Equating the (e^^^ + e*^^^) components results in 



1 . 



4/1I/2 



(F.5) 



and the e''^^ component gives 

^2g-2k 



hi 



Ah 



2b3<^' - 3/163$' - 4/1635' - 2/i6'3 



+ KNce--^^'h2 (a^ - 2a6 + l) . (F.6) 



This leaves 62 and 63 to be determined. Substituting these results into 
(F.4), we find that the (e''^^ + e'"^^) component of H-^ = dB2 reduces to the 
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equation of motion (B.9) for b. The only remaining equation is then the 
gpSv component. This is a first order differential equation in 62 and 63, 

= 8he^3+^% + 2 (a^ - l) + e^^^+'^^h' Ua^ - l) e^%3 + 4e^% 



(F.7) 



Solving for 62 we obtain 



SV/l 



(a^ - 1) e4f+2'^/i'63 



4aa' + (4^' + - 4^' - 



3 ' 



(F. 



which does not appear to be very useful. Instead, we can use the fact that 
we want Qpage, D3 = (see eq. 6.12). We therefore impose that the 6^"^^^^ 
component of — B2 A -F3 vanishes. The resulting equation is algebraic in 
62 and 63, and results in 



-2h 



4/1I/2 



e2»/ii (1 - a2) b: 



K 3* 

3 e 2 



7V2(a_6)5' + 4e2(9+'^)$' 



(F.9) 



Together with the above results for 61^4^5 this completes (6.2). 



It remains to check that this 62 is also compatible with the requirement 



that dB2 = H3. Substituting into (F.7) we find that 63 cancels, giving 
= 4e^(9+/^) 1 2/1 



A 2 



2g'^' + 2h'<^' + <l>" + 2 {<!>') 



2g'h' - 2h'h' - h 



"} 



a^e^s - 2a%e^3 + 2(0 - b)b"e^^3+h) + 4^^ _ 5)5'e2(g+h)^/ 



+ 2abe^3 _ 2 {p'f e^i9+h) _ ^ig _ 



Ah 



(F.IO) 
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This is solved by the equations of motion (B.7, B.9) for <I> and b. 

To determine the effect of the undetermined function 63, we can look at 
the difference AB2 = B2 — (-62)63=0) which we find to be of the form 

A52 = Fi{p) sin 61 d9 Ad(f + F2{p) sin9 dO A d(f + Fs{p) cos 6* dp A dip 

+ F4^{p)cose dp A dip + F^ip) dp A dip, (F.ll) 

where the Fi depend on g, h, 63 and their derivatives. If we set this equal 
to 



d 



l3i{p) cos 9 dp + /32{p) cos 9 dip + /33{p) dip 



(F.12) 



we can solve for the giving 
1 



AB2 



-d 



2g+*/2 ^ f^f^^ (pQg 9 d^ + COS 9 dp + dip) 



(F.13) 



G Appendix: Seiberg-like duality 



In section 6.1.6 we discuss how the operation known as Seiberg duality in 
the KS cascade acts for our non-SUSY solution. In order to do so, we find it 
instructive to compare to two different cases: the KS case and the baryonic 
branch case. These are summarized here. 



G.l The KS case 

We follow here the treatment in j20j. specified in the case of no flavors 
{Nf = 0). The NS potential B2 is given by, 

^2 = ^[/<?l A52 + %3 A(74] (G.l) 

where the definition of (71, 54 can be found in [2QJ. When specialized to 
the cycle 

^2 = [9 = 9,ip = 2TT -ip,ip = iPo] (G.2) 

we obtain that 

B2\^2 = ^[{f + k) + {k - f)cosipQ\sm9d9 Ad^ (G.3) 
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from which one finds 



2 ' 



(G.4) 



On the other hand, as computed in [20j, we can see that the Maxwell charge 
of D3 branes is 



QMax,D3 = —[f-{f-k)F] 

vr 



We see that under the change 



vr 



k-^k 



IT 



the D3-Maxwell charge changes by 

QMax,D3 QMax,D3 " ^c, ^0 — ^ ^0 — 1- 



(G.5) 



(G.6) 



(G.7) 



these transformations, are equivalent to changing the NS potential with a 
large gauge transformation 



B2 + -[gi A 52 + ff3 A gi] 



(G. 



which when evaluated on the cycle S2, produces the changes in eq.(G.7). 
We move now to analyze the baryonic branch (SUSY) case. 



G.2 Baryonic branch case 

In this case the NS potential is 



B2 



Ke 



/lV2 



[eP^ - cos a{e^^ + e^^) - sina(e^2 j_ g<pi)] 9) 



Evaluating this on the S2 we get 

,2$ 



vr 



+ /) + (A;-/)cosVo 



k + f 
k-f 



Ke 



2<l> 



N, 
Ke 



=2g 



COS a(^(a^ + 1) - e^'') + sin Qae''+^' 



c '- 

2'!' r g2g 



COS a a + sin ae^~^^ 



(G.IO) 
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Using the explicit expressions, we have 

2<I> 2^ 

~k = -'^QcotHp), / = -^Qtanh(p) (G.ll) 
The Maxwell charge for D3 branes can be written as, 

QMa.,D3 = -e2^+2h+2$ci>' (G.12) 



vr 

and using the BPS equation for ^' we have 

vr 



N 

QMax,D3 = — 2f + {k- f)F 



(G.13) 



where F = {1 — b). So, once again, we obtain that under a large gauge 
transformation, 

bo bo — 1, QMax,D3 QMax,D3 — ^c- (G.14) 
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